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We present a measurement of the cross section for W-boson production in association with jets 
in pp coUisions at \/s = 1.96 TeV. The analysis uses a data sample corresponding to an integrated 
luminosity of 320 pb~^ collected with the CDF II detector. W bosons axe identified in their electron 
decay channel and jets are reconstructed using a cone algorithm. For each W+ > n-jet sample 
(n = 1 — 4) we measure the differential cross section da-(j>p W-\- > n-jet)/di5^*''~^*' x B{W 
ev) with respect to the transverse energy Et of the n"'-highest Et jet above 20 GeV, and the total 
cross section a{pp W-\- > n-jet; EJ^*'"'^"* > 25 GeV) x B{W -> ei^), for a restricted W ^ eu 
decay phase space. The cross sections, corrected for all detector effects, can be directly compared 
to particle level W + jet(s) predictions. We present here comparisons to leading order and next-to- 
leading order predictions. 

PACS numbers: 14.70.Fm, 13.87.Ce, 12.38.Qk, 13.85.Ni 
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Final states containing a vector boson V {V ~ W, Z) 
and multiple jets (y-l-jet(s)) are a key signal channel for 
important standard model (SM) processes such as tt or 
single top production, as well as a search channel for the 
Higgs boson and for physics beyond the SM. The produc- 
tion of y-|-jet(s) via quantum chromodynamics (QCD) 
presents a very large background to these processes. The 
ability to describe it accurately is therefore crucial, as 
well as being a stringent test of the power of perturbative 
QCD predictions. Consequently, a precise measurement 
of the cross section for QCD V -I- jet(s) production is an 
important component of the hadron collider experimen- 
tal program. In this paper, we report a measurement [ij 
of the differential cross sections for direct W ev+ > n- 
jet production as a function of the transverse energy 
S^*"* y of the n*^-leading jet (the highest Et jet for 
W+ > 1-jet, the second highest Et jet for W+ > 2- 
jet, etc.), for n = 1 - 4 and E^^''* > 20 GeV. We also 
provide the total cross section dn — a{W ev+ > n- 
jet; j^1^*^~3^* > 25 GeV) for n = 1 — 4. In order to min- 
imize the dependence of the measurement on the mod- 
eling of the W boson production and decay kinematics, 
we quote cross sections defined within a limited W decay 
phase space: E^i^ > 20 GeV, iTy^''^] < 1.1, E!f > 30 GeV, 
and > 20 GeV/c^ [2]. The range of S^"^* extends 
up to 350 GeV in the W + 1-jet sample, a significant 
increase in the measured phase space compared to pre- 
vious y-l-jet(s) measurements Furthermore, the 
differential spectra presented here are for the first time 
corrected for all detector effects and represent absolute 
particle level cross sections Q free, within systematic un- 
certainties, of any experimental bias. As such, they pro- 
vide a benchmark which can be directly used for back- 
ground estimates and for the validation and tuning of 
QCD phenomenological models. At the end of this pa- 
per we present, as examples, comparisons of our results 
with some of the available predictions. It is important 
to note that the cross section is not corrected for effects 
resulting from the interaction between the proton and 
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anti- proton remnants (the "underlying event"). Such a 
correction would introduce into the measurement a de- 
pendence on theoretical models of the underlying event. 

This analysis uses 320 ± 18 pb~^ of data collected us- 
ing the CDF II [0] detector during the Tevatron Run II 
period. The CDF II detector is a general-purpose de- 
tector designed to study pp collisions at the Fermilab 
Tevatron. Inside a 1.4 T solenoidal magnetic field, a 
large open-cell drift chamber and an eight-layer silicon 
system provide precise charged-particle tracking informa- 
tion. Outside the solenoid electromagnetic and hadronic 
sampling calorimeters surround the tracking volume, al- 
lowing for the measurement of particle energies over the 
range \ri\ < 3.6. Finely segmented detectors located at 
electromagnetic shower maximum are used for electron 
identification. Forward gas Cerenkov detectors measure 
the fraction of bunch crossings that result in an inelas- 
tic pp collision and thereby determine the instantaneous 
luminosity delivered to the experiment. 

The selection oi W ^ ev events proceeds as follows. 
An online trigger system selects events containing an elec- 
tromagnetic calorimeter cluster with Et > 18 GeV asso- 
ciated with a high pT track. Offline, electron candidates 
are required to pass standard identification cuts Q and 
to have E^'^ > 20 GeV. The sample is enriched with 
events containing a neutrino by requiring that the miss- 
ing transverse energy I^t 0, corrected for the jet energy 
scale (see below) and the potential presence of muons [s!] , 
satisfies I^t > 30 GeV. To reduce background contami- 
nation further, the W transverse mass [3| rn^ is required 
to satisfy > 20 GeV/c^. In addition, Z -> e+e" 
events are rejected by a veto algorithm |7i] . 

The jets in each W ^ ev event are reconstructed 
using the jetclu cone algorithm Q with cone radius 
R — ^ A(/)^ + /Sjf- = 0.4. Starting from seed locations 
corresponding to calorimeter towers with Et > 1 GeV, 
all nearby towers with Et > 0.1 GeV are used to search 
for stable cones. To resolve ambiguities with overlapping 
cones, cones sharing an energy fraction greater than 0.75 
are merged into a single jet; otherwise the shared towers 
are assigned to the closest jet. We apply a jet energy scale 
(JES) correction [l3| such that the measured Ei^* is on 
average equal to the summed Et of the particles within 
the jet cone that are the result of the pp ^ W + X in- 
teraction. Jets are required to have E^* > 20 GeV and 
hi < 2.0. 

To ensure negligible overlap of electron and jet en- 
ergy deposits, events are rejected if any jet lies within 
AR = 0.52 of the W decay electron [1]. Using this jet 
definition, we divide the inclusive W ev candidate 
events into W+ > n-jet samples, and in each sample 
form the E^* spectrum of the n*'' i?T-ordered jet. 

The processes which contribute background events to 
our W candidate sample can be divided into two cate- 
gories: "leptonic" and "multi-jet". The leptonic back- 
ground contains real electrons and/or neutrinos from bo- 
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son decay and includes W tv, Z — > e+e~, WW, W^, 
and top pair production. The multi-jet background arises 
from QCD interaction events in which one or more jets 
are incorrectly reconstructed in the detector as electrons 
and have mis-measured energy, resulting in large event 
I^T- Background estimation proceeds as follows. For 
each W+ > rt-jet sample, a background-enriched event 
sample is constructed by removing the I^t > 30 GeV 
requirement in the W selection. Multi-jet, leptonic, and 
signal I^T histograms are then fit to the data in the range 
[0,100] GeV. The measured I^t spectrum, and the result 
of the fit, are shown in Fig. [l]for the W+ > 1-jet sam- 
ple. For this fit the leptonic background and signal pro- 
cesses are modeled by applying the W event selection 
minus the I^t > 30 GeV requirement to detector simu- 
lated Monte Carlo event samples of these processes. The 
multi-jet background is modeled using an event sample 
selected from the same 320 pb^^ analysis dataset by re- 
quiring that at least two of the electron identification 
criteria fail. Kinematic cuts are unchanged, resulting in 
a background-dominated sample that accurately reflects 
the kinematic distributions of the multi-jet background 
events in the signal sample. It is necessary to correct 
this multi-jet sample for ^ 5% contamination from sig- 
nal and leptonic background events, estimated by apply- 
ing the multi-jet selection criteria to the Monte Carlo 
simulations of these processes. In the fit, only the nor- 
malizations of the multi-jet and signal I^t histograms 
arc allowed to float. The normalization of the W ^ tv 
and Z e+e^ histograms relative to the signal is fixed 
by the well-established relationships between these cross 
sections The normalizations of the WW, W^, and 
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top pair production histograms are determined using the 
recently measured cross sections for these processes [l^ . 
Once the fit is performed, the background fractions in 
each W+ > n-jet sample are obtained by integrating 
the respective histograms, with their fitted normaliza- 
tion, above the I^t cut of 30 GeV. These fractions are 
then used to normalize the jet Et distributions of each 
background model relative to the candidates to give the 
background correction as a function of Ei!^^ . This method 
offers a more accurate description of the kinematics of 
the multi-jet background when compared with previous 
approaches 3J. Fig. [T] demonstrates a successful model- 
ing of the I^T spectrum. Similarly good agreement was 
found in the electron Et and mjf distributions across all 
jet multiplicities. 

Rarely, reconstructed jets may originate from separate 
pp interactions in the same bunch crossing. To account 
for this effect, corrections to the measured cross sections 
are computed by multiplying the number of overlapping 
pp interactions, estimated using the primary vertex mul- 
tiplicity in the signal sample, with the rate of jet produc- 
tion in "minimum-bias" events selected independently of 
activity in the central detector. The corrections are less 
than 2% in the lowest Ei^* bins and decrease rapidly with 
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FIG. 1: Tlie results of fitting the signal and background 
I^T template distributions to the data in the W+ > 1-jet sam- 
ple before the final I^t cut is applied. 



increasing Ei^^ . 

The total background fraction increases with increas- 
ing jet multiplicity and transverse energy. At low Ei^^ 
it is 10% (40%) in the 1-jet (4-jet) sample, rising to 90% 



at the highest Ei^* for all jet multiplicities. Multi-jet 
events contribute ~ 70% of the overall background in 
the 1-jet sample. At high jet multiplicities and high E^* 
the contribution from top pair production becomes in- 
creasingly important, climbing to 50% (80%) of the total 
background in the 2-jet (3, 4-jet) sample. The systematic 
uncertainty on the background estimate is 15% at low 
-Ey*^* independent of the jet multiplicity, rising to 50% 
(20%) at the highest E^''* in the 1-jet (4-jet) sample. At 
low jet multiplicities, this is dominated by the limited 
statistics of the multi-jet background sample. At high 
jet multiplicities, the 12% uncertainty on the measured 
top pair production cross section dominates the total sys- 
tematic. 

We have used fully simulated signal Monte Carlo sam- 
ples to correct the event yield for the efficiency of the 
W —I- ev selection criteria. Samples for each jet mul- 
tiplicity n were obtained using the ALPGEN vl.3 (isj 
event generator for the W + n-paiton final state, and the 
PYTHIA v6.3 IJ] Monte Carlo program for the parton 
shower and hadronization. PYTHIA includes an underly- 
ing event model, hereafter referred to as tune a, which 
has been tuned to describe Tevatron data [l3|- The ef- 
ficiency correction factor is defined as the ratio of two 
subsets of the generated signal events: in the numera- 
tor, the number of reconstructed events which pass the 
W eiy selection criteria, and in the denominator, the 
number of generator-level events which pass the electron, 
neutrino and transverse-mass cuts corresponding to our 
cross section definition. This is found to be (60 ±3)%, in- 
dependent of event jet multiplicity and kinematics within 
the quoted uncertainty. By comparing Z — s- e+e^ mea- 
sured and simulated event samples, we confirm that elec- 
tron identification efficiencies are well reproduced by our 
Monte Carlo across all measured phase space to within 
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5%. This uncertainty also covers the observed variation 
in efficiency obtained by changing the number of final 
state partons in the ALPGEN matrix element (ME) calcu- 
lation and the parton showering program from PYTHIA to 
HERWIG v6.5 [la]- Additionally, the W candidate event 
yield must be corrected to account for the efficiency of the 
online trigger to accept high Et electrons. This is inde- 
pendent of jet kinematics and found to be 96.2±0.6% 

A further correction to the event yield is required to 
form particle level W — > et^-f > n-jet cross sections as a 
function of Ei!^* . This correction factor accounts for the 
effect of calorimeter jet energy resolution on the mea- 
sured cross section and is determined as follows. Using 
the ALPGEN-PYTHIA simulated signal samples, two cross 
sections in each Ei!j^^ bin are determined: one defined by 
clustering generator-level particles into jets, the other by 
clustering after detector reconstruction, using the same 
calorimeter level jet definition (including JES correction) 
as the one used in the data. The correction factor for 
each bin is then defined as the ratio of the particle to 
calorimeter level cross section in that bin. To avoid de- 
pendence of such a correction on the assumed particle 
level E^* distribution, an iterative procedure is used to 
reweight the events at the particle level until the parti- 
cle level E^* distribution agrees with the corresponding 
data-unfolded distribution to within the systematic un- 
certainties of the measurement. The correction factors 
vary between 0.95 and 1.2 over the measured range of 
E't'\ 

The total systematic uncertainty on the cross section 
introduced by the jet energy measurement ranges be- 
tween 5% and 20%, increasing with increasing Ei!^* . This 
is dominated by the approximately 3% [lOj] uncertainty 
on the JES correction. The effect of this uncertainty on 
the cross section is estimated by applying this variation to 
an ALPGEN-PYTHIA simulated signal sample reweighted 
to match the data. The sensitivity of the measurement 
to jet energy resolution uncertainties, estimated by vary- 
ing the calorimeter resolution in the simulation, is much 
smaller by comparison. 

The measured differential cross sections da{W 
ev+ > n-jet)/(i£'^*''~-''^* are listed in Table HI For each 
inclusive jet multiplicity sample the cross sections are 
given with respect to the Et of the n*''-leading jet in 
the sample, i?^*''"-''^'. The quoted statistical uncertain- 
ties are on the event yield in each bin. The system- 
atic uncertainties are the sum in quadrature of the ef- 
fects introduced by the uncertainty on the background 
estimation, acceptance correction, and jet energy mea- 
surement. A 5.8% uncertainty in the integrated lumi- 
nosity is not included, since this uncertainty is com- 
pletely correlated between different bins. In 
summary, the total systematic uncertainty on the mea- 
sured cross sections is < 20% at low E^r^* increas- 
ing to 50% - 80% at high Ei^ 



dominated by the uncertainty on the jet energy scale, 
whereas at high Ei^* it is dominated by the background 
uncertainty. We also provide the total cross section 
cr„ = a{W -> ev+ > n-jet; i;^*''"^'^* > 25 GeV) for 
n = 1-4: 0-1 53.5±0.6(stat.)±4.6(syst.)±3.1(lum.)pb; 
CT2 = 6.8 ± 0.2(stat.)±1.0(syst.)±0.4(lum.) pb; erg = 
0.84± 0.10(stat.)±0.21(syst.)±0.05(lum.) pb; and 0-4 = 
0.074 ± 0.039(stat.)±0.035(syst.)±0.004(lum.) pb. The 



choice of Er 



ith—jet 



> 25 GeV is made in order to pro- 



,jet 



At low -Ey*^* this is 



vide a benchmark measurement that is less sensitive 
to the impact of the underlying event, largest at low 
Ei^* (see below). We include for completeness the to- 
tal inclusive pp ^ W x B{W ev) cross section 
for the restricted W ^ ev decay phase space: cto = 
798± 2(stat.)±40(syst.)±46(lum.) pb. 

We proceed to compare the measured cross sections 
to some of the available theoretical predictions. Leading 
order (LO) perturbative QCD calculations exist for the 
matrix element of ^ -I- n partons, with n < 6. They are 
included in Monte Carlo event generators [13, [13, 
where the initial and final state partons are evolved 
through a perturbative parton shower (PS) and eventu- 
ally hadronized. Additionally, the generator may include 
a model of the underlying event. In this LO plus PS 
approach, ambiguities may arise as a result of the hard 
emission of gluon radiation during the parton shower evo- 
lution. For example, a y -f n-parton event may be recon- 
structed, after the shower, with a jet multiplicity rij ^ n, 
and the question naturally arises as to whether the event 
should or should not be counted in the estimate of the 
rij-jet cross section. This problem has been studied ex- 
tensively in the literature, leading to the development of 
three merging algorithms, usually known as CKKW [l9| . 
Lonnblad's [20|, and MLM [21] . A merging algorithm 
ensures that a given configuration in the multi-jet phase 
space enters into the calculation once and only once. 

We present here the first comparisons of -I- jet (s) im- 
plementations of the CKKW and MLM schemes to data. 
We use an implementation of the CKKW scheme here- 
after referred to as SMPR '^3\ . Details of the generation 
parameters and systematic uncertainties are given in 22 1 
for the SMPR model and in [il for the MLM model. The 
former uses madgraph v4 [l7| for the ME generation, 
PYTHIA v6.3 for the PS, and CTEQ6L1 ^ parton distri- 
bution functions (PDFs), the latter uses ALPGEN v2.12, 
HERWIG v6.5, and CTEQ5L [25*] PDFs respectively. Fol- 
lowing generation, the SMPR and MLM predictions are 
formed by clustering the final state particles into jets us- 
ing the JETCLU algorithm. The uncertainties on these 
predictions cover variation of the renormalization scale 
by a factor 0.5-2. This dominates the overall uncertainty 
in the absolute rates [2^] . 

In addition, we present comparisons of our data with 
next-to-leading order (NLO) predictions for -I- 1 and 
W + 2 jets obtained with the MCFM [26] program. 
These were generated using the CTEQ6.1M PDFs 
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FIG. 2: Top: the ratio of data to theory for the total cross 
sections o-„ = cr(W^ ev+ > n-jet; > 25 GeV) as a 

function of the jet multiphcity n. Bottom: an/(Jn~-i for data, 
MLM, SMPR and MCFM. Inner (outer) error bars denote the 
statistical (total) uncertainties on the measured cross sections. 



and a renormalization and factorization scale 



+ (Pt')^- We define an uncertainty due to the 

jet 

Pt ' 



choice of /i by generating with a lower scale, /i 

Addition- 



and a higher scale, /i = 2 * 



'"W 



ally, the variation due to the uncertainty on the PDFs 
has been computed using the Hessian method 2j| . This 
PDF uncertainty is also broadly applicable to the SMPR 
and MLM predictions. 

In the case of the NLO predictions, the final states 
are not evolved through a parton shower nor hadronized. 
Jets are reconstructed with a cone algorithm R = 0.4, 
such that two partons are merged if they are within 
1.3x i? of each other and within R of the resulting jet cen- 
troid [22,] • This is still considered to be sufficient to give 
a reasonable description of the perturbative structure of 
the jet ^28^. However, before comparing with data, the 
non-perturbative effects of hadronization and the under- 
lying event have to be considered. We have estimated, 
using PYTHIA TUNE A, the impact of these two effects. 
The effect of the underlying event is to increase the cross 
section with respect to the parton level, whilst the effect 
of hadronization is to decrease it. The magnitude of both 
effects decreases asymptotically with increasing E^r^*" . Be- 
low 50 GeV the hadronization effect dominates, leading 
to an overall decrease of the cross section with respect to 
the parton level that is within 10%. At higher E^r^*" , the 
correction is driven by the underlying event leading to 
an increase of at most 5%. A detailed study of this cor- 
rection is outside the scope of this paper, and we do not 
apply any such corrections to the MCFM predictions. 

The upper plot of Fig. [2] shows, as a function of the 
jet multiplicity n, the ratio of data to theory for the to- 
tal cross sections ct„ — (j{W ev+ > n-jet; e^^~-''^^ > 
25 GeV). The lower plot shows the ratio (T„/(t„_i. In 
Fig.[3]the ratios of the measured differential cross sections 
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FIG. 3: Ratio of the measured cross sections da{W —> 
eu+ > n-jet)/d£;y*''"^'=' to the MLM, SMPR and MCFM 
predictions for n = 1 (top), n — 2 (middle) and n = 3 (bot- 
tom). MCFM predictions are not available for n = 3. Inner 
(outer) error bars denote the statistical (total) uncertainties 
on the measured cross sections. 



da{W ev+ > ?i-jet)/d£^-^*''~-''^* to the predictions are 
shown for 71 = 1 — 3. The difference observed in Fig. [2] 
between the measured cross sections and SMPR or MLM 
predictions reflects the LO nature of these calculations. 
All the predictions show good agreement with the data 
in the cross section ratios (T„/cr„_i. Fig.[3]shows that the 
variation in the W + n-jet cross section as a function of 
Ei^* is better reproduced by the SMPR prediction than 
by the MLM. A possible explanation is the absence of a 
tuned underlying event model in the HERWIG component 
of the MLM prediction. We observe good agreement be- 
tween the MCFM predictions and data in both total and 
differential cross section comparisons. 

In summary, we have used 320 pb~^ of CDF II data to 
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TABLE I: The measured cross section da{W ei'+ > n- 
jet)/d£^"'--"=* for each Sj"^--''=*-bin (n = 1-4), with statis- 
tical (first) and systematic (second) uncertainties. An overall 
5.8% uncertainty on the integrated luminosity is not included. 





da/dE^'*-''* (pb/GeV) 


20- 


25 


4.46 ± 0.07 ±0.29 


25- 


30 


2.80 ± 0.06 ±0.21 


30- 


35 


1.92 ± 0.05 ±0.17 


35- 


40 


1.31 ± 0.04 ±0.14 


40- 


50 


0.839 ± 0.023 ±0.093 


50- 


60 


0.498 ± 0.018 ±0.063 


60- 


75 


0.259 ± 0.011 ±0.038 


75- 


90 


0.158 ± 0.008 ±0.024 


90- 


110 


0.056 ± 0.005 ±0.011 


110 


-150 


0.0225 ± 0.0023 ±0.0044 


150 


-195 


0.0035 ± 0.0011 ±0.0023 


195 


-350 


0.00039 ± 0.00022 ±0.00038 


— tt; 

H/rp 


(GeV) 


da/dE^ ^ (pb/GcV) 


20- 


25 


0.874 ± 0.033 ±0.097 


25- 


30 


0.483 ± 0.025 ±0.066 


30- 


35 


0.286 ± 0.020 ±0.045 


35- 


40 


0.190 ± 0.017 ±0.034 


40- 


50 


0.095 ± 0.009 ±0.019 


50- 


60 


0.057 ± 0.007 ±0.012 


60- 


75 


U.UZOy It U.UU'iO ztU.UUDD 


75- 


95 


n m n7 -u n nnoo -un nnoc 
U.UiU( =t U.UUzz ztU.UUzo 


95- 


190 


0.00059 ± 0.00039 ±0.00049 


H/rp 


"-'^^ (GeV) 


da/dE':^''"^'^* (pb/GeV) 


20- 


25 


0.184 ± 0.016 ±0.036 


25- 


30 


0.087 ± 0.012 ±0.024 


30- 


35 


0.037 ± 0.008 ±0.013 


35- 


45 


0.020 ± 0.006 ±0.011 


45- 


80 


0.0015 ± 0.0012 ±0.0013 




da/dE^'"-^'-' (pb/GeV) 


20- 


25 


0.0422 ± 0.0087 ±0.0079 


25- 


35 


0.0074 ± 0.0039 ±0.0036 



measure the differential cross section da{pp W+ > n- 
]et) / dEl^^~-''^* X B{W ev) as a function of the trans- 
verse energy of the n*^-leading jet, and the total cross 
section cr„ = ^(pp ^ W+ > n-jet; E'^^''^^"^ > 25 GeV) x 
B{W eu), for n = 1 — 4 in a restricted W ^ ev decay 
phase space. The cross sections, corrected for all detec- 
tor effects, can be directly compared to the particle level 
predictions of W + jet(s) Monte Carlo generators. 
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